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^ ■ Abstract. Recently, the diffusion-limited cluster aggregation (DLCA) model 

■ was restudied as a real-world example of showing discontinuous percolation 
' transitions (PTs). Because a larger cluster is less mobile in Brownian motion, 

it comes into contact with other clusters less frequently. Thus, the formation of a 
giant cluster is suppressed in the DLCA process. All clusters grow continuously 
with respect to time, but the largest cluster grows drastically with respect to 

■ the number of cluster merging events. Here, we study the discontinuous PT 

occurring in the DLCA model in more general dimensions such as two, three, and 

' four dimensions. PTs are also studied for a generalized velocity, which scales with 

cluster size sasvs oc . For Brownian motion of hard spheres in three dimensions, 
the mean relative speed scales as s~^^'^ and the collision rate aVs scales as ~ s^^^. 

■ We find numerically that the PT type changes from discontinuous to continuous 
(-H ' as T] crosses over a tricritical point rjc ~ 1.2 (in two dimensions), jjc ~ 0.8 (in 

■ three dimensions), and rjc ~ 0.4 (in four dimensions). We illustrate the root of 

this crossover behavior from the perspective of the heterogeneity of cluster-size 
distribution. Finally, we study the reaction-limited cluster aggregation (RLCA) 
model in the Brownian process, in which cluster merging takes place with finite 
probability r. We find that the PTs in two and three dimensions arc discontinuous 
even for small r such as r = 10"'^, but are continuous in four dimensions. 
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1. Introduction 

The notion of percolation is widely used to explain the formation of a macroscopic 
spanning cluster in diverse systems [Tj. In percolation, as the control parameter, i.e., 
the density of occupied nodes (in site percolation) or bonds (in bond percolation) , is 
increased, a macroscopic spanning cluster emerges at the percolation threshold. This 
behavior is referred to as the percolation transition (PT), which is conventionally 
continuous. More generally, the term PT is used for the emergence of a macroscopic 
cluster in growing networks, and is occasionally referred to as PT. In the random graph 
model introduced by Erdos and Renyi (ER) [5], the control parameter is the number 
of bonds (links), and a macroscopic giant cluster is found to emerge at a critical point. 
That is, PT occurs at a finite percolation threshold. Recently, Achlioptas et al. [3] 
proposed a modified ER model, in which the growth of the largest cluster is suppressed 
drrring the dynamical evolution. In this model, the macroscopic giant cluster emerges 
at a delayed transition point and the transition occurs in a rather explosive manner. 
Thus, the PT in this model is claimed to be discontinuous. Following the proposal of 
this explosive percolation model, many studies have been performed on discontinuous 
PT; however, whether such explosive percolation transitions are indeed discontinuous 
in the thermodynamic limit is still a matter of debate and sensitive to detailed dynamic 
rules [lEllilllllinilinilllllllllllll]. Nonetheless, the introduction of such 
explosive percolation models has led to intensive studies of discontinuous PT. In this 
circumstance, we wonder whether such discontinuous PTs indeed can be observed in 
real- world systems. 

In our previous work |15j . we studied the diffusion- limited cluster aggregation 
(DLCA) model following Brownian motion in two dimensions for simplicity, as an 
example of a real-world system showing discontinuous PT. Here, the PT means the 
formation of a giant component, as conventionally used in the evolution of random 
graphs, instead of the formation of a spanning cluster, as conventionally used in regular 
lattice. We monitored the PT as a function of the number of cluster aggregations. 
Because real-world systems can be three dimensional, here we extend our previous 
study to three- and four-dimensional cases. Moreover, we study the PT for the case of 
cluster velocity being in the general form Vs oc s"^ and find that the PT type changes 
from discontinuous to continuous as ry increases. In the last part of this paper, we 
extend our study to the reaction-limited cluster aggregation (RLCA) model in which 
clusters diffuse following Brownian motion, and when two clusters come into contact 
with each other, they merge with a certain probability r and remain separate with the 
remaining probability 1 — r. We find that the discontinuous PT behavior can also be 
observed in this RLCA model in two and three dimensions but that the PT remains 
continuous in four dimensions. 

2. DifFusion-limited cluster aggregation model 

The DLCA model was introduced by Meakin et al. [16] and Kolb et al. [l^ . Initially, 
A'^ particles are placed randomly in a d-dimensional lattice space with linear size L. 
The density of the particles is fixed as p = N/L'^, whereas the system size L is 
controllable. Simulations start from N monoparticles. When an s-sized cluster moves 
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in the Brownian process, its velocity is given as ~ s~^/^, and the collision rate 
per cluster becomes av^ ^ s^/^ in three dimensions when the cluster is regarded as 
a hard sphere [TH]- This originates from the fact that a Brownian particle with mass 
m has mean velocity v = and velocity fluctuations ksT/m when the particle 
is in thermal equilibrium with temperature T, where the overbar means ensemble 
average over thermal fluctuations and fcs is the Boltzmann constant. Thus, we obtain 

's/v^ ^ l/^/rn, which leads to Vs ~ s~^/^ when cluster mass is regarded as being 
linearly proportional to cluster size. Accordingly, the choice of a cluster of size s with 
probability s-^/^/ i^g^^^^ ^ s'^/^. 

To implement this velocity in simulations, we perform a simulation in the following 
steps: Initially, all particles are single, (i) An s-sized cluster is picked up with a 
probability proportional to s^^/^, and it is moved to one of the nearest-neighbor 
positions. All particles in a mobile cluster move together with the cluster shape 
unchanged. After this move, when two distinct clusters come into contact, those 
clusters merge with probability one, forming a larger cluster, (ii) Time is advanced by 
St = ^/Cl2s NsS~^^^), whenever the cluster moves irrespective of whether a contact 
occurs, but the control parameter p is advanced by 1 /N only when the cluster is placed 
next to another cluster and the two clusters merge. When all particles merge into a 
single cluster, the dynamics ends. 

We presume that cluster aggregations take place irreversibly. Then the number of 
cluster merging events during the whole process is — 1. For example, if two immobile 
clusters are merged by one mobile cluster, and thus the three clusters become one, 
then the number of merging events is counted as two, and p is advanced by 2/N. 
Thus, the number of merging events corresponds to the number of inter-cluster edges 
connected. When a created cluster contains loop structure, then the number of cluster 
merging events is not same as the number of occupied bonds in bond percolations. 

The parameter p represents the number of cluster merging events per total particle 
number, and corresponds to the number of links connecting two distinct clusters per 
network size in the random graph model. The variable p turns out to differ from time t 
in a nontrivial way. The order parameter of the PT is the giant cluster size per system 
size N, denoted as Gn = G. To examine the PT. Gn is measured as a function of p, 
which becomes nonzero beyond the transition point pc, where PT is discontinuous if 
G^ipt) ^ Gn{Pc) > ^iiid continuous if it is zero. 



2.1. Brownian motion 

We begin by recalling our previous work in two dimensions. The giant cluster size 
increases monotonically as a function of time t; however, it increases drastically when 
monitored as a function of the variable p as shown in Figsma) and 1(b), respectively. 
This indicates that the difference originates from the nonlinear relationship between t 
and p shown in Figlljc). In particular, the time interval between two successive cluster 
merging events becomes long when p approaches one, because few clusters remain and 
they hardly ever contact each other. Thus, such a nontrivial relationship between t 
and p arises. 

To verify the discontinuity of the order parameter, we use the finite-size scaling 
approach, which is different from the conventional one used for the continuous PT 
|19j . In this approach, a particular point Pd{N) was introduced as a triggering 

p,iN)^p^{N)~{^^\j"G^{p.), (1) 
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Figure 1. (Color online) (a)-(f) Simulation results of the DLCA model in two 
dimensions with veloeity exponent r) = —0.5 for Brownian motion, (a) Plot of 
G vs t. G grows monotonically from t = with inereasing t. (b) G grows 
discontinuously near p 1 with respect to p. (c) Relation between t and p. t 
grows rapidly as p approaches p = 1, by which a discontinuous PT occurs with 
respect to p. (d)-(f) Finite-size scaling analysis for the discontinuous PT of G, xii 
and X2 using Eqs. I I2I4I I. Data of different-size systems collapse well onto a single 
curve, which implies that a discontinuous PT occurs. We use the system sizes 
= 6, 10, 14, 18, 22, and 26 with density p = 0.05 for numerical simulations, 
(g) Numerical estimations of k^/C' at pd for rj = — 0.5(0), ri = 0(o), rj = 0.4(A), 
and rj = 0.8(V). Simulations are performed in systems with TV = 8000 and 
L = 400. Slopes of the guidelines are -0.23 ± 0.02, 0.32 ± 0.04, 0.62 ± 0.01, 
and 0.88 ± 0.01 from above, (h) Numerical estimations of ki/C with the same 
symbols as in (g) used for each rj. Here = 8000 and L = 400 are used for 
simulations. Slopes of the guidelines are 0.35 ± 0.04, 0.32 ± 0.04, 0.2 ± 0.01, and 
0.05 ± 0.01 from above, (i) Plot of G vs p of different size systems with r; = 1.5 in 
two dimensions. A giant component emerges continuously near p ^ 0.45. Data of 
AT/IO^ = 8,32, 72, 128, and 200 with p = 0.05 are shown in this plot, (j) Finite- 
size scaling analysis for continuous PT for the data used in plot (i). pc(oo) = 0.43, 
l/v = 0.29, and ji/v = 0.5 are estimated. 
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where Pc{N) is the point at which the slope of G^lp) becomes maximum. It is found 
that dG jsi {p) / dp\p^ increases in a power-law manner '-^ A^^/" with l/F « 0.86 ± 0.02. 
Then, since the giant component size grows as GNiPf) — GM{Pd) ^ ^(1): where pf 
is the final step of cluster aggregation, during the interval pj — Pd ^ 0{N~'^/'^), the 
transition is indeed discontinuous. The above behavior is also checked by using the 
scaling ansatz for the discontinuous PT, 

G^(p) (xiV-/^/Vo((p-p,(iV))iVi/^), (2) 

where P/v = for the discontinuous transition. This scaling form differs from that 
conventionally used for continuous transitions, which is written as Eq. (jH]) shown later. 
Thus a discontinuous PT can be confirmed by checking whether the data of Gm{p) 
versus /o(?3) with p = [p — pd)N^^'^ for different system sizes collapse onto a single 
curve or not. Indeed, we find that the data from different system sizes collapse onto 
a single curve in two dimensions when using the value l/v ~ 0.86 in Fig. [Ijd). 

We also studied the susceptibility, defined as xi(p) = Y^'s ^'^''^s{p) / Y^'^ sns{p), 
where the prime represents the exclusion of the giant component in summation. This 
function can be represented in the scaling form, 

Xi{p)^N'<-'^f,{{p-pd)N^/^). (3) 

It was found that the data from different system sizes collapsed well onto a single 
curve with the exponent value 71/F = 0.95 in Fig. [Ije). 

We also attempt a scaling analysis for another quantity of the susceptibility X2 
defined as X2 = ^ \/ {{G{p) — (G(p)))^). This quantity is the standard deviation of 
G{p) for a given p. We can check that X2 also collapses well onto a single curve, 

X2(p)«^^^/V2(b-Pd)iV'/^). (4) 

with 72/^7 = 1 in Fig. [ijf). This result suggests that the PT is indeed discontinuous. 

Similar analyses are carried out in three and four dimensions. In three dimensions, 
we obtain results similar to those of two dimensions but with different exponent values, 
i.e., 1/V = 0.70 ± 0.02, 71/77 = 0.90 ± 0.02, and 72/17 ^ 1.00 ± 0.02. The data from 
different system sizes collapse well onto the scaling functions, which are shown in 
Fig. [2] In four dimensions, we obtain similar results but with different exponent 
values, i.e., 1/17 = 0.60 ± 0.02, 7i/l7= 0.92 ± 0.01, and 72/17= 1.00 ±0.01. The data 
from different system sizes collapse well onto the scaling functions as shown in Fig. [3l 
These results verify the discontinuity of PT in three and four dimensions. 

We also investigate the discontinuous PT in a different approach via the 
Smoluchowski equation, which describes the dynamics of cluster aggregations. In 
particular, we introduce an asymmetric Smoluchowski equation in which the collision 
kernel is different depending on whether cluster is mobile [TS] as follows: 

dn^ ^i^j n^kg risks . . 

where Us = Ns/N is the concentration of s-sized clusters, and k'^jG' and kjjG are 
the collision kernels of immobile and mobile clusters, respectively. C" = k'sUs and 
C = 'Y^^kgUs are the normalization factors. The first term on the right-hand side 
of Eq. ([5]) represents the aggregation of a mobile cluster of size i and an immobile 
cluster of size j with s = i + j. The second term represents a mobile cluster of 
size s merging with an immobile cluster of any size including the largest cluster, in 
which kjHj ~ G(jj) is used. The third term represents an immobile cluster of 
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Figure 2. (Color online) Finite-size scaling analysis for discontinuous PT of 
the DLCA model in the Brownian process (r; = —0.5) in three dimensions 
having linear sizes L/10 = 25,30,35,40,45, and 50 with p = 0.004. (a) Data 
collapse of G by using Eq. l[2} for different size systems, (b) The estimation 
of the exponent 1/F defined in the relation dG(p) / dp\p^ ~ N^^" . We obtain 
1/F = 0.70 ± 0.02. (c) Plot of xi for different system sizes in a scaling form, (d) 
The estimation of the exponent 71 /F defined in the relation xi{Pcl) ~ TV'^i'''. 
We obtain 71 /F = 0.90 ± 0.02. (e) Data collapse of X2 of different system sizes 
in a scaling form, (f) The estimation of the exponent 72 /F defined in the relation 
X2{Pc2) ~ N^i^/'^. We obtain 72/F = 1.00 ± 0.02. In particular, we can find 
that the peak height of the standard deviation of the largest cluster increases 
proportionally to the system size, which is a feature of the first-order transition. 
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Figure 3. (Color online) Finite-size scaling analysis for discontinuous PT of the 
DLCA model in Brownian motion {r) = —0.5) in four dimensions. We use the 
scaling forms of Eq. l(2} for G in (a), Eq. Q for xi in (c), and Eq. |4]for X2 in 
(e). Here, the systems having linear size L/10 = 8, 10, 12, 14, 16, 18, and 20 with 
p = 0.0004 are used for simulations. These data are well collapsed in each scaling 
form, which confirms that the PT is indeed a discontinuous PT. (b), (d), and (f) 
Numerical estimations of the exponents used in each scaling form. In particular, 
we can find 72/57 = 1, which is a feature of the first-order transition. 



size s merging with a mobile cluster of any size including the largest size, in which 
^'j^j ~ ^'{p) used. The summation runs only for finite clusters. Once an 
infinite-sized cluster is formed, the dynamics is terminated. We also do not need to 
consider finite clusters and infinite cluster separately as in sol-gel transitions [50] . The 
collision kernel is determined by intuitive argument as the perimeter of clusters and 
thus kyC ^ i^+^-^/df for mobile and ki/C ^ i^-'^l'^s for immobile clusters, where 
dj is the fractal dimension of clusters [21] and r\ = —0.5 for the Brownian case. 

In two dimensions, using df « 1.4, it is estimated that k[ ^ i^^ "^ and ki ~ 
which are in agreement with the numerical estimations of k'^ ^ and ki ^ 

shown in Figs, mg) and 1(h), respectively. In three dimensions, dj ~ 1.8 [52], and 



thus fc' 



-0.1 



and ki 



MA 



are expected, which are again in agreement with the 
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Figure 4. (Color online) Plots of (a) k[/C' and (b) ki/C vs i at of the DLCA 
model in three dimensions. The system ol N = 4000 and L = 100 is used, (a) 
The estimated slopes are -0.08 ± 0.03 {r) = -0.5) (□), 0.43 ± 0.02 (-q = 0.0) (o), 
0.85 ±0.01 (»? = 0.4) (A), and 1.35 ±0.01 {-q = 1.0) (v). (b) The estimated slopes 
are 0.43 ± 0.03 (□), 0.45 ± 0.03 (o), 0.34 ± 0.03 (A), and 0.20 ± 0.03 (v). 
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Figure 5. (Color online) Plots of (a) k[/C' and (b) ki/C vs i at p^j of the DLCA 
model in four dimensions. Simulations arc carried out in the system oi N = 5184 
and L = 60. (a) The estimated slopes are 0.04 ±0.02 (?} = -0.5) (□), 0.64 ±0.02 
(r? = 0.0) (o), 0.97±0.04 (»? = 0.4) (A), and 1.33±0.02 (r? = 0.8) (v) from above, 
(b) The estimated slopes are 0.64 ± 0.03(n), 0.64 ± 0.02{o), 0.57 ± 0.02(A), and 
0.38 ± 0.04(V). 
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Figure 6. (Color online) Finite-size scaling analysis for the discontinuous PT 
based on the numerical simulation results of the Smoluchowski equation. Collision 
kernels obtained from the data in Fig. [l] Fig. [4] and Fig. [5] for the Brownian 
case are used in two, three, and four dimensions, respectively. The system 
sizes we used are N/IO'^ = 16,32,64,128,256,512,1024, and 2048. (a)-(c) In 
two dimensions, dG(pc)/dp behaves as ~ N. By using this result, data collapse 
behaviors are obtained for (a) the giant cluster size, and the susceptibilities (b) xi 
and (c) X2- (d)-(f) Similar analyses are carried out for three dimensions, in which 
dG{pc) / dp ~ N. (g)-(i) Similar analyses are carried out for four dimensions, in 
which dG{pc)/dp ~ 



measured values k'^ ~ z^"-^ and ki ^ as shown in Fig. |4l In four dimensions, 
df ~ 2.0 [22], and thus fc- w const and ki ~ are expected, which are again in 
agreement with the measured values fc^ ^ i'^''^ and ki ^ as shown in Fig. [5l 

Next, we investigate the growth of the giant component size G by simulating the 
Smoluchowski equation numerically with the collision kernels we measured. Starting 
from N monomers initially, numerical simulations are carried out as a function of p for 
different system sizes. We plot the giant cluster size G{p) and the susceptibilities xi 
and X2 as a function of p in scaling forms, and we find that the data of different system 
sizes collapse well onto a single curve with the critical exponents previous obtained as 
shown in FiglS] 

2.2. Generalization of velocity scaling 

In this section, we study the PT of the DLCA model. Here, the scaling of the collision 
rate is generalized for computational simplicity by shifting the scaling with cluster 
size entirely into the scaling of the cluster velocity: Vs ^ s'^. It is then necessary to 
know whether the scaling exponent of the cluster velocity rj can be positively valued. 
Consider the motion of a small solid sphere introduced into a non-uniform electric field 
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Figure 7. (Color online) Snapshot of clusters in a three-dimensional lattice at 
p = 0.997 with the system size TV = 4000 and L = 100. There are eleven clusters 
for each snapshot. We find that cluster sizes are more heterogeneous as rj grows. 
This means that the rate of growth of the largest cluster increases as ri increases. 

in air. This arrangement is readily realized in a corona discharge |23] . for example, 
between a grounded hollow cylinder and a thin wire along the cylinder's axis, when 
the wire is charged to a negative high voltage with respect to the cylinder. The corona 
current sustains drifting negative ions towards the cylinder walls. The sphere attracts 
the ions onto its surface by the image charge effect until the Coulomb repulsion by 
the accumulated ions prohibits it. The sphere is accelerated by the local electric field 
while its motion is resisted by the Stokes drag [53], reaching a terminal velocity that 
scales as the radius of the sphere. In 3-D the mean velocity scales as Vs ^ s^^^ , and 
the collision rate scales as avs ~ s for solid spheres. For fractal spheres, as in the 
DLCA model, r] can become positive. 

To implement the effect of this velocity in simulations, we select an s-sized cluster 
with a probability proportional to s"^ , allow the cluster to diffuse to a nearest neighbor, 
and make time pass by 5t = '^/i^^NsS^). If the cluster comes into contact with 
another cluster, the variable p is advanced by 1/-/V, regardless of the cluster size 
s. As rj increases in positive region, the velocity of large clusters becomes large, so 
that they have higher probability of merging with another cluster. Thus, the growth 
rate of larger clusters is higher. This behavior was originally observed by Meakin et 
al. [12]. They argued that the exponent of cluster size distribution increases as the 
velocity exponent increases. In this case, the PT is continuous, because the giant 
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Figure 8. (Color online) Plots of the cluster size distribution ns{pt) ~ s~'^ for 
several ri around the tricritical point r)c. The slopes of the guidelines are (a) —1.8, 
—2.0, and —2.3 in two dimensions, (b) —1.7, —2.0, and —2.4 in three dimensions, 
and (c) —1.7, —2.0, and —2.4 in four dimensions. The system sizes we used are 
(a) Af = 8, 000, L = 400, (b) N = 13, 500, L = 150, and (c) = 40, 000, L = 100. 
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Figure 9. (Color online) (a) Plot of the giant cluster size G vs p of the 
DLCA model with velocity exponent T] = 1.0 > »)c ~ 0.8 in three dimensions. 
Data are obtained by Monte Carlo simulations, (b) Numerical simulation of the 
Smoluchowski equation. In (b), we use the collision kernel fc^/C ~ i^-^^ and 
ki/C ~ i"'^*^. Inset in (a): Data collapse of the data plotted in (a). Linear sizes 
of the systems are L/10 = 10, 15, 20, 25, 30, 35, and 40 and the density of particles 
is p = 0.004. Data of different system sizes collapse well onto a single curve 
predicted theoretically for a continuous transition in which we use pc(oo) = 0.56, 
1/v = 0.2, and j3/v = 0.6. The behavior of the data collapse indicates that the 
transition is indeed continuous. Inset in (b): Data collapse of the data plotted 
in (b). We use system sizes N/W* = 8,32,128,512, and 2048 and substitute 
Pc(oo) = 0.66, 1/v = 0.14, and fijv = 0.55 for data collapse. Again the data- 
collapse behavior indicates a continuous PT. 
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cluster grows continuously. In contrast, when -q is negative, large-sized clusters are 
suppressed in growth, and their number is reduced. Instead, medium-sized clusters 
become abundant. As p increases, such medium-sized clusters merge suddenly and 
create a giant cluster. Thus, the PT is discontinuous. Fig. [7] shows the snapshots of 
clusters for different values of r/ just before the percolation threshold. From these 
properties, one can guess that the transition type changes from discontinuous to 
continuous as rj increases across a certain value ?7c- 

To determine the tricritical point ?7c, we start from rj = —0.5 and observe the 
change of the cluster size distribution by increasing rj. In our previous study |25] . it 
was shown that discontinuous (continuous) PT occurs when ns(p^) ~ s~'^ satisfies 
T < 2 (t > 2). We use this result to determine the tricritical point rjc- Fig. [8] shows 
the cluster size distribution for several values of ry near jy^. We estimate r/^ ~ 1.3 in 
two dimensions [Fig. [SJa)], w 0.8 in three dimensions [Fig. Mjo)] and w 0.4 in four 
dimensions [Fig. [SJc)] based on numerical data. 

To confirm the continuity of PT in the region rj > r/c, we perform finite-size 
analysis for G{p). The behavior of G{p) with 77 = 1.5 > rjc ~ 1.3 in two dimensions is 
plotted in Fig.[lji). As we can observe in Fig.[lji), the crossing point of G between two 
different sizes decreases as the system size increases, which means that the transition is 
continuous. But this tendency cannot be seen clearly in this figure. Thus we attempt 
a finite-size scaling analysis for continuous transition. In Fig. [TJj), the G{p) of the 
different system sizes used in Fig. [TJi) is collapsed onto a single curve in the scaling 
form. 

Gn{p) oc N-^Z-'foiip - p,{oo))N^'''). (6) 

A similar analysis is carried out in three dimensions, which is shown in Fig. |9|a). In 
this analysis, we use the system of 77 = 1.0 > rjc ~ 0.8. In the inset, we use the 
scaling form given by Eq. ([6]). To verify the continuous transition in an alternative 
way, we use the Smoluchowski equation. Similar to what we found in the previous 
subsection, we obtain k[/C' ~ i^'^^ and fci/C « at pd in Fig. S] Fig. IHtb) shows 
the simulation result of the Smoluchowski equation for various system sizes. In the 
inset, we use the scaling form given by Eq. (|6]) to verify continuity and find that data 
are well collapsed on a single curve. These results confirm that the transition is indeed 
continuous in the region 77 > 77c « 0.8 in three dimensions. 

A similar analysis used for three dimensions is also applied for four dimensions, 
which is shown in Fig. [TOl In Fig. rTOT a). the G{p) of various system sizes are plotted 
and these data are well collapsed onto a single curve if we use the previous scaling 
form given by Eq. ([6|), which is shown in the inset. Second, we take k[ ^ i^'^*^ and 
ki ^ i^'^^ from Fig. [5] to simulate the Smoluchowski equation in Fig. rTOT b'). In the 
inset, the data are well collapsed in the previous scaling form Eq. IB] These results 
again confirm the continuous transition in the region rj > rjc ~ 0.4 in four dimensions. 

3. Reaction-limited cluster aggregation model 

Here, we perform similar studies for the reaction-limited cluster aggregation (RLCA) 
model in the Brownian process in two, three, and four dimensions. In this model, two 
adjoining clusters merge irreversibly with probability r, but with the remaining with 
probability 1 — r, they can move independently. As r goes to 0, a cluster can penetrate 
inside the area between branches of another cluster, becoming trapped and irreversibly 
stuck within it. As a result, the resulting cluster becomes less ramified, and its fractal 
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Figure 10. (Color online) (a) Simulation result of DLCA with velocity exponent 
ri = 0.8 > T]c ~ 0.4 in four dimensions, (b) Numerical simulation of the 
Smoluchowski equation using the collision kernel k'^/C' ~ i^-^^ and ki/C ~ i"'*". 
Inset in (a); Data collapse of the data used in (a). We perform simulations with 
linear sizes L/10 = 6,8,10,12, and 16 with particle density p = 0.0004. Data 
of different system sizes collapse well onto a single curve predicted theoretically 
for a continuous transition in which we use pc(oo) = 0.56, 1/v = 0.22, and 
P/u = 0.43. The behavior of the data collapse indicates that the transition is 
indeed continuous. Inset in (b): Data collapse of the data used in (b). Sizes 
of the systems are AT/IO" = 8,32,128,512 and 2048. We use pc(oo) = 0.64, 
1/v = 0.25, and fi/v = 0.3 for data collapse. Again the data-collapse behavior 
indicates a continuous PT. 
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(P-P<1~)W°'^ (P-P,(~))W"^^ Ip-P^.(~)IA*' 



Figure 11. (Color online) (a)-(d) Plot of the RLCA model in the Brownian 
process in two dimensions, (a) Plot of and ki/C vs i at p^. N = 2000 

and L = 100 are used in the simulations. The two collision kernels are estimated 
as k'JC ~ j-0.l2±0.06 and ki/C ~ jO.SOio.OS (b)-(d) Plots of finite-size scaling 
behaviors of (b) G, (c) Xli s-id (d) X2- simulations of the Smoluchowski equation 
are performed using N/10'^ = 16,32,64,128,256,512,1024, and 2048. Data 
are well collapsed onto the scaling curves for the discontinuous transition with 
dG{pc)/dp ~ Af in two dimensions, which indicates that the PT of the RLCA 
model in the Brownian process is discontinuous. Plots (e)-(h) correspond to plots 
(a)-(d) but in three dimensions. For (e), N = 4000 and L = 100 are used. The 
collision kernels are estimated as k'JC' ~ j-0.04±0.03 g^^d ki/C ~ i0.70±o.04^ 
In (f)-(h), the system sizes are Af/10* = 16,32,64,128,256,512,1024 and 2048. 
Data collapse is well behaved using the formulas for the discontinuous transition 
with dG(pc) /dp ~ A'"'''^^. Plots (i)-(l) correspond to the plots (a)-(d) but in four 
dimensions. For (i), A^ = 5184 and L = 60 are used. The collision kernels 
are estimated as k'-/C' ~ j0.29±0.0l ^nd ki/C ~ j0.90±0.00^ {3)-{l), the 

system sizes are A/10-* = 16, 32, 64, 128, 256, 512, 1024 and 2048. Data are well 
collapsed on the scaling formulas for the continuous transition. It is estimated 
that Pc(oo) = 0.95, 1/u = 0.21, and i3/u = 0.05. 



dimension is increased |26| . Here we use r = 10~^, in which the dynamics of cluster 
aggregations of the RLCA observed is different from that of the DLCA. 

To obtain the giant cluster size G, we measure the collision kernels k[/C' and ki/C 
at pd and simulate the Smoluchowski equation using these collision kernels. Monte 
Carlo simulations of the RLCA model take extremely long computation times for us to 
understand the finite-size scaling behavior of the PT. Thus, we measure the collision 
kernels in two, three, and four dimensions, and then investigate the finite-size scaling 
behavior of numerical data of the Smoluchowski equations. The measured collision 
kernels are shown in Fig. 1111 The collision kernel may be written in a power-law form, 
k[kj - Then the exponents (w', w) are estimated as (-0.12 ±0.06, 0.50 ±0, 03), 

(-0.04 ± 0.03,0.70 ± 0.04) and (0.29 ± 0.10,0.90 ± 0.00) in two, three and four 
dimensions, respectively. It is noteworthy that in the conventional Smoluchowski 
equation, clusters are immobile, and thus the collision kernel is symmetric as kikj. 
If kikj ^ {ij)'', then for ( > 0.5, the percolation transition is continuous. In this 
case, the cluster-size distribution follows a power law at the critical point as Us ^ 
where r > 2 |20j . However, for the asymmetric case above, the criterion for continuous 
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transitions is not known specifically to our knowledge. We find that the cluster size 
distribution for the asymmetric Smoluchowski equation for the RLCA model with the 
numerically estimated kernels at the transition point follows a power law, ~ g-o-so 
in two dimensions, Us ^ s~^''^^ in three dimensions up to finite-size cutoffs, and 
Us ^ 5-2.05 jjj^ fp^j. dimensions. Thus, the percolation transition for four dimensions 
can be expected to be continuous. We also remark that the collision kernel for the 
RLCA does not agree well with the one obtained from the formulas k[ ^ f^-'^/d-f+'n 
and ki ^ . If we use df = 1.7, df = 2.0, and df = 2.4 for two, three and four 

dimensions, respectively, and 77 = —0.5 [l?]- This difference is because merging of 
two clusters does not occur at their perimeters in the RLCA process. Using these 
obtained collision kernels, we perform numerical simulations of the Smoluchowski 
equation, and we find that the giant cluster size G and the susceptibilities xi ^-nd 
X2 behave following the scaling functions for the discontinuous PT in Eqs.(I2]|3]) in two 
and three dimensions. However, in four dimensions, the obtained data do not collapse 
onto the scaling functions of the discontinuous transition, but rather collapse onto 
Eq. (|6]), valid for continuous transitions. This different behavior is caused by the large 
exponent values of the collision kernels. 

4. Summary 

In this paper, we extended the previous study of discontinuous PT of the diffusion- 
limited cluster aggregation (DLCA) model to three and four dimensions. We showed 
that the discontinuous PT also occurs even in three and four dimensions for Brownian 
motion. In this case, the discontinuous PT is caused by the natural suppression effect 
of Brownian motion to the growth of large clusters. Moreover, we studied PT for 
the DLCA model with general velocity Vs ^ for various values of 77, where s is 
the cluster size. As 77 increases, the suppression effect becomes weak, so that there 
exists a tricritical point rjc, across which the PT type changes from discontinuous 
to continuous. Finally, we briefly studied the PT for the reaction-limited cluster 
aggregation (RLCA) model in Brownian motion in two, three and four dimensions. By 
simulating the Smoluchowski equation with the obtained collision kernels, we find that 
the PT is discontinuous in two and three dimensions but continuous in four dimensions. 
In this work, Pc ^ 1 for the cases of discontinuous transitions, otherwise Pc < 1 in 
the cases of continuous transitions. Conclusively, we expect that the discontinuous 
PT can be observed in many modified DLCA models owing to the suppression effect 
of Brownian motion. 
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